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A COHERENT TIMING SOLUTION FOR THE NEARBY, THERMALLY EMITTING ISOLATED NEUTRON 
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ABSTRACT 

We present a phase-coherent timing solution for RX J0420. 0—5022, the coolest (kT ~ 45 eV) and 
fastest-spinning (P = 3.45 s) of the seven so-called isolated neutron stars (INSs). Using 14 observations 
with the XMM-Newton spacecraft in 2010-2011, we were able to measure a spin-down rate v = 
(-2.3 ± 0.2) x 10" 15 Hzs _1 (P = (2.8 ± 0.3) x 10" 14 s s" 1 ), from which we infer a dipolar magnetic 
field of 1.0 x 10 13 G. With reasonable confidence we were able to extend the timing solution back 
to archival XMM-Newton from 2002 and 2003, giving the same solution but with considerably more 
precision. This gives RX J0420. 0—5022 the lowest dipole magnetic field of the INSs. Our spectroscopy 
does not confirm the broad absorption feature at 0.3 keV hinted at in earlier observations, although 
difficulties in background subtraction near that energy make conclusions difficult. With this, all 6 of 
the INSs that have confirmed periodicities now have constrained spin-downs from coherent solutions. 
The evidence that the INSs are qualitatively different from rotation-powered pulsars now appears 
robust, with the likely conclusion that their characteristic ages are systematically older than their true 
ages, because their fields have decayed. The field decay probably also causes them to be systematically 
hotter than pulsars of the same (true) age. 

Subject headings: stars: individual (RX J0420. 0—5022) — stars: neutron — X-rays: stars 



1. INTRODUCTION 



Th e so-c alled iso l ated n eutron stars (INSs; see iHaberU 
120071 and iKaplanl 120081 for reviews) are a group of 
seven nearby (< 1 kpc) neutron stars with low (~ 
10 32 e rgs" 1 ) X-ray luminosities and lo ng (3-10 s) spin pe- 
riods (Kaplan fc van Kerkwi ik 2009bj). They are unique 
in that the X-ray emission likely comes from a large 
fraction of the neutron stars' surfaces, and is not in- 
fluenced by accretion (as in the case of X-ray binaries) 
or non-thermal magnetospheric emission (as in the case 
of rotation -powered pulsars); the INSs are also radio- 
quiet (e.g.. iKondratiev et al.ll2009l ). The INSs then have 
the potential to help understand neutron star radii and 
cooling via measurements of their emission areas and 
luminosities, but this is made difficult by our inability 
to realistically model the X-ray, ultra- violet, and op- 
tical emission from these objects (e.g., IHo et al.l 120071 : 
IKanlan et alll201lh . 

Rece n tly it was proposed dKaplan fc van Kerkwiik 
20095 iPons. Miralles. fc Geppertl 120091: iPopov et al ' 



the broad absorption features seen at en ergies of 0.2- 
0.75 keV in the spectra of alm ost all INSs (|Haberll 120071: 
van Kerk wiik fc Kaplan|[2007l ). also requires knowledge 
of the magnetic fields, as at these field strengths the as- 
sumed transition energies are field-dependent. There- 
fore, we have undertaken systematic measurements of 
the dipole magnetic fields for the INSs through phase- 
coherent X-ray timin g using the Chandra and XMM - 
Newton spacecraft ([K aplan & van Kcrkwii kl I2005al[g ; 



van Kerkwiik fc Kaplan! 120081: i Kaplan fc van Kerkwiil 
2009allbl hereafter IKvK05al : IKvK05bt IvKKOSt IKvKOff 



KyK09b). 



2010T ) that the current temperatures and magnetic fields 



of the INSs reflect non-thermal, coupled evolution, where 
the magnetic field has decayed in strength, heating the 
neutron-star surface. Testing this hypothesis is an impor- 
tant step to using the INSs to constrain the overall cool- 
ing history of neutron stars, and th rough it probe their 
inner structure and composition ( |Yakovlev fc Pethickl 
2004). To constrain evolutionary models, measure- 
ments of the current magnetic fields and temperatures 
of the INSs are required. Similarly, to understand 
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Here we measure the spin-down of the INS 
RX J0420.0-5022 (hereafter RX J0420), the last INS 
with a confirmed period measurement but without a 
timing solution. RX J0420 was identified as a p ossi- 
ble neutron star bv lHaberl. Pietsch. fc Motchl (|1999h . al- 
though it had originally been associated with a nearby 
galaxy. Followup ROSAT observations showed a very 
soft thermal spectrum and yielded an improved posi- 
tion, both of which led to its classification as a neutron 
star. While initially a 22.7-s period was suggested, ob- 
servations with XMM found a period of 3.45 s in stead, 
the shortest among all INSs (jHaberl et al.l [2001 here- 
after IH+Q4h : the same observations also showed that 
RX J0420 was the coolest INS, with kT ~ 45 eV. 

2. OBSERVATIONS & ANALYSIS 

We observed RX J0420 fourteen times with XMM 
(| Jansen et all 120011) in 2010 and 2011, and focus here 
on the data taken with the European Photon Imaging 
Camera (EPIC) with pn and MOS detectors, all used in 
small window mode with thin filters (Tabled]). We re- 
processed our observations with SAS version 11.0.0 and 
calibration files current as of 2011 May 25. We also re- 
processed the pn data from lH+04 which are taken with 
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TABLE 1 

Log of Observations and Times of Arrival 



2004 



2006 



2008 



2010 



2012 







Exp. a / bg a 


TOA b 


>, 


Rev. 


Date 


(ks) Counts 3 (%) 


(MJD TDB) 





560 


9009 Hop ^0 
j_^e^ ou 


20 047 


4 201 


10.0 


561 


9009 Dpr ^1 


20 048 


4 593 


11.7 


570 




20 547 


4 647 


11.6 


664 


2003 Jul 25 


20,036 


4,384 


11.8 


1887 


2010 Mar 30 


7,472 


2,062 


41.2 


1890 


2010 Apr 04 


9,072 


2,679 


45.8 


1892 


2010 Apr 09 


7,772 


2,049 


35.4 


1913 


2010 May 21 


5,471 


1,462 


32.1 


1948 


2010 Jul 29 


6,472 


1,626 


37.2 


1975 


2010 Sep 21 


9,872 


2,372 


38.2 


1981 


2010 Oct 02 


11,672 


2,980 


34.8 


1981 


2010 Oct 03 


12,972 


3,064 


38.6 


1981 


2010 Oct 04 


16,871 


4,465 


40.2 


1983 


2010 Oct 06 


10,471 


2,586 


36.9 


2008 


2010 Nov 26 


5,471 


1,339 


35.2 


2032 


2011 Jan 13 


15,471 


4,223 


44.5 


2071 


2011 Mar 31 


7,018 


1,577 


38.2 


2076 


2011 Apr 11 


5,471 


1,248 


32.2 



52638 
52640 
52658 
52846 
55285 
55290 
55295 
55337 
55406 
55460 
55472 
55472 
55473 
55476 
55526 
55575 
55651 
55662 



2855519(12) 

0466236(12) 

8319656(8) 

0226435(10) 

5413317(17) 

8432112(23) 

4037368(11) 

2761247(20) 

6382597(30) 

4233059(20) 

0349984(22) 

8839787(11) 

3194015(13) 

0218532(19) 

4316848(19) 

0260821(17) 

9068660(24) 

3338203(19) 
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Fig. 1. — Phase residuals for RX J0420. In the top panel, we show 
the residuals relative to a linear model (y = 0). The line shows 
the best-fit quadratic solution. Residuals relative to the quadratic 
solution are shown in the bottom panel. Inset: pulse profile for 
RX J0420, based on a combination of all of the 2010-2011 data 
folded according to the ephemeris in Tables [2] and corrected for the 
background. The blue curve is the best-fit sinusoidal profile. 



Note. — All observations used the small window mode and thin 
filter for both EPIC-pn and EPIC-MOS1/2, except for Revs. 560, 
561, 570, 664, in which the full window mode was used (which 
meant that only the EPIC-pn data were suitable for timing). 
a The exposure time, number of counts, and estimated fraction of 
events due to background f^ s given here are for EPIC-pn only. 
b The TOA is defined as the time of maximum light of the fun- 
damental closest to the middle of each observation computed from 
the combined EPIC-pn and EPIC-MOS1/2 datasets, and is given 
with 1-0" uncertainties. 

the same filter, but with the full window mode instead 
(we did not use their full-frame MOS data, since these 
do not resolve the pulsations). We used epchain and 
emcliain. and selected source events from a circular re- 
gion of 37'.'5 radius. For the pn, we selected energies 
between 130 and 800 eV, where we set our lower energy 
cutoff slightly below the default of 150 eV to increase 
the net number of counts from our very soft target (by 
about 25%; for even lower thresholds, the instrumental 
background increases too rapidly), and the upper cutoff 
at the energy at which the source becomes undetectable 
(thus minimizing the effects of flares, which dominate 
the background at higher energies). Following standard 
practice, we included only one and two-pixel (single and 
double patterns 0-4) events with no warning flags for pn, 
and single, double, and triple events (patterns 0-12) with 
the default flag mask for MOS1/2. We barycentered the 
event tim es usin g the Chandra X-ray Observatory posi- 
tion from|H±04j: a = 04 h 20 m 01^95 and S = -50°22'48'.'l 
(J2000). We extracted background lightcurves for pn 
from similarly sized regions offset from the source, but 
at the same RAWY coordinate, as recommended by the 
SAS User Guide0 For MOS1/2, the small- window mode 
does not permit such large background areas, but we 
used several smaller areas to compensate. 

2.1. Timing Analysis 

Our timi ng analy sis largely follows the procedure de- 
scribed in [KvK05al As a starting place, we first de- 
termined the frequency that maximized the power in 
a Z\ periodogram for the EPIC-pn data from the 

3 See http://xmm.esac . esa. int/external/xmm_user_support/do 



longest observation in Rev. 1981. We then expanded 
the periodogram to include data from all observations 
in Revs. 1981 and 1983, finding a best-fit frequency of 
v = 0.289 6033±0.0000003Hz, consistent with that found 
by IH+04I for the earlier data. In contrast to some of the 
other INSs, there is no evidence for higher harmonics in 
the periodogram: the Z\ power is 33.8, while Z\ = 34.7 
and -Zf = 35.2, both of which are consistent with the ad- 
ditional power of 1 expected for noise. We also checked 
to see if the true period was in fact 6.9 s (closer to that 
of the other INSs), but the pulse shape, hardness ratio, 
and median energy were all consistent with no variation 
between the first and second halves of the pulse (with 16 
bins the lightcurve variation between the first and sec- 
ond halves has \ 2 — 10.6/8, while the hardness ratio 
variation has \ 2 — 4.1/8). 

Using the above frequency, we determined the times-of- 
arrival (TO As; see Table Q} for the combined EPIC data 
from each observati on by fitting the binned lightcurves 
(following KvK05b) to a single sinusoid, appropriate 
given the results of the periodograms (Fig.[T]); the best- fit 
sinusoid to the composite, background corrected pn data 
had a semi-amplitude of 15 ± 2%, where the uncertainty 
includes an estimate for the variation in the background 
correction over different background regions. This is 
a litt le high er than the semi-amplitude of about 10% 
from IH+04L but differences in background-subtraction 
and energy selection could account for the difference (our 
best-fit to the 2002-2003 data has a semi-amplitude of 
13.2 ± 1.1%). The x 2 f° r the fit to the composite profile 
was good, 14.4 for 13 dof. 

Using our TOAs, we were able to identify a reasonably 
unambiguous coherent timing solution. This was possible 
as we restricted solutions to have < 9 x 10~ 13 Hzs -1 
or BcH n < 2 x 10 14 G (based on the incoherent limits 
set by lHabe"rl 120071 ) . Among those, the solution pre- 
sented in Table [2] was the best, yielding \ 2 = 19-9 for 
11 degrees of freedom, with alternatives at \ 2 = 25.1 
(B di p = 1.8 x 10 14 G), x 2 = 37.3 {B dip = 7 x 10 13 G), 
and x 2 — 38.5 (z> > 0). Of these, all but the first can 
be excluded on statistical grounds. The first comes from 
c ™ ei ftfl\ffl?;e%%aiM®-/gf ^L^f^k^khe cycle count between the 
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TABLE 2 

Measured and Derived Timing Parameters for 
RX J0420. 0-5022 



Quantity 




Value 






2010-2011 




2002-2011 



Dates (MJD) . . 


55,286-55,662 


52,638-55,662 


t (MJD) 


55,430.6001387(6) 


55430.6001387(6) 


v (Hz) 


0.2896029061(12) 


0.2896029058(10) 


v (lO" 15 Hzs" 1 ) 


-2.3(2) 


-2.314(8) 


TOA mis (s) ... 


0.3 


0.2 


X 2 /DOF 


19.9/11 


24.8/15 


P(s) 


3.453004024(14) 


3.453004027(12) 


P (lO-^ss" 1 ). 


2.8(3) 


2.759(10) 


fchar (Myr) 


2.0 


2.0 


B dip (10 13 G) . . 


1.0 


1.0 


E (10 31 ergs -1 ) 


2.7 


2.7 



Note. — Quantities in parentheses are the formal 1-cr uncer- 
tainties on the last digit. t c ^ t = P/2P is the characteristic age, 
assuming an initial spin period Pq -C P and a constant magnetic 

field; -Bnin = 3.2 X 10 19 V PP G is the magnetic field inferred as- 



3 dip 

suming spin-down by dipole radiation; E - 
the spin-down luminosity. 



3.9 X 10 4( W 



erg s 



densely sampled Rev. 1981-1983 group and the next clos- 
est obs e rvatio n, Rev. 1975, and is close to the limit from 
lHaberll (|2007h : we will return to this alias shortly. We 
were able to identify the same solut ion usin g a single co- 
herent Z\ (y, v) periodogram (as in|vKK0l|. Spin-down 
is well-detected, at ^ 10 cr. The reduced \ 2 1S somewhat 
high, but even adjusting our uncertainties to allow for 
an reduced \ 2 of 1 will still give an 8cr detection of spin- 
down. The implied magnetic hel d is well w ithin the range 
of other detections for the INSs (|KvK09b[ ). 

We can confirm and improve our solution by extrapo- 
lating it back to the older data from 2002-2003. The 
time difference is roughly 2600 days and our v uncer- 
tainty gives a formal cycle-count uncertainty of ±5 cycles, 
but by trying multiple solutions, we find that only a sin- 
gle cycle count difference leads to a solution that fits all 
four earlier TOAs. Trying this generally, iteratively ex- 
ploring all cycle-count ambiguities between all data sets, 
we find a single best-fit solution that agrees with the 
best-fit solution using only the new data (Table [5] and 
Figure [IJ. With \ 2 = 24.8 for 15 degrees-of-freedom, it 
still has a slightly high reduced \ 2 ■ However, the alter- 
nate solutions either have nearly the same implied spin- 
down rate but differ slightly in the cycle counts between 
the 2002-2003 and 2010-2011 observations (and the low- 
est of those has \ 2 = 33.0), or are dramatically different 
but have significantly worse \ 2 (x 2 = 45.9 and i> > 0). 
The alternate, £?dip = 1-8 X 10 14 G solution from above 
does not extrapolate well to the older data, allowing us to 
exclude it. Overall, the older data thus allow us to con- 
fidently select the correct solution to the new data, and 
select a reasonably secure overall solution (A% 2 = 8.1). 

2.2. Spectroscopic Analysis 

We examined all EPIC-pn spectra of RX J0420. (A 
full spectral analysis, including the EPIC-MOS and RGS 
data and a phase-resolved analysis, is in progress.) We 
used the same source and background extraction regions 
as for the timing analysis, created appropriate response 




Fig. 2. — Blackbody fits to the merged EPIC-pn data of 
RX J0420, from 2002-2003 (blue circles) and 2010-2011 (red 
squares). The model is the thick black line. The middle panel 
shows the backgrounds used, and highlights the background struc- 
ture seen at low energies in the 2010-2011 data. The bottom panel 
shows the residuals. 

files, and binned the spectral files such that the number 
of source plus background counts was at least 25 and the 
bin width was at least 30 eV (so that there are roughly 
2 bins per EPIC-pn resolution element). 

We first compared the raw EPIC-pn spectra of all of 
the observations against each other. This did not include 
any response files or calibration corrections, but even so 
the binned pn spectra were generally consistent with each 
other, implying no spectral change (Fig. [2]). There are 
small deviations at low energies (0.2-0.4 keV) and we will 
return to these below, but the 2010-2011 data did not 
show any appreciable variability. 

We fit the pn data using sherpa (Ref sdal et al.| [2009). 
To aid in fitting we merged the event and response files 
into t wo gro ups: one for 2002-2003 (full-frame data, also 
fit bv lH+041 ) and one for 2010-2011 (small-window data). 
While not perfect, we found that an absorbed blackbody 
provided a reasonable fit, with Ah < 1 x 10 18 cm~ 2 , 
kToo = 47.6 ± 0.3 eV, and = 12.8 ± 0.3kmkpc _1 
(formal 1-cr uncertainti es; y 2 = 73.0 for 31 dof), rea- 
sonably consistent with IH+04I . Unlike IH+04L we do not 
find evidence for a spectral feature, perhaps because of 
changes in the response files and calibration since ear- 
lier fits. However, there are indications that our fit is 
not completely reliable. First, the best-fit value of the 
absorption is (although it is covariant with the black- 
body temperature). Secon d, there are some residuals 
near 0.33 keV, where lH+041 found evidence for a spectral 
line. However, we could not find a consistent fit to both 
sets of data. Third, stronger differences are seen at en- 
ergies of 0.2-0.3 keV (see Figure [2]). Some of these are 
too narrow to come from astrophysical sources, and in- 
stead likely reflect problems in background subtraction 
(the low-energy background for the small-window data 
in particular can be significant and has substantial en- 
ergy structural). Hence, it is difficult to interpret any 
residual structure there with respect to a blackbody. 

3. DISCUSSION & CONCLUSIONS 

We have determined a reliable, statistically significant 
coherent spin-down solution for RX J0420. With this, 

4 See http : //www . star . le . ac . uk/$\sim$amr30/BG/BGTable . html 
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only RX J1605. 3+3249, which as yet has only a ten- 
tative detection of a periodicity, lacks a coherent solu- 
tion (although in the cases of RX J2143. 0+0654 and 
RX J0806.4— 4123 spin-down was not well measured, and 
further observations are in progress). While the over- 
all results of our timing program have been disc ussed at 
length in previous pape rs (over all energetics in [yKK08; 
spectral implications in lKvK09t evolutionary models in 
lKvK09bl ). here we touch on some of the aspects that 
make RX J0420 unique and compare it to the other ob- 
jects in its class. 

First, while the timing properties (dipole magnetic 
field, E, characteristic age) of RX J0420 place it well 
within the INSs, RX J0420 has the shortest period by 
more than a factor of 2. In the context of the magneto- 
thermal evolution model, this could be a consequence 
of a lower initial magnetic field, and thus less dramatic 
early spin-down. RX J0420 also has the lowest current 
field, although it is not clear whether there is a good 
correlation between current magnetic field and period: 
RX J0806.4-4123 and RX J2143.0+0654 both have long 
periods but relatively weak magnetic fields. The mag- 
netic field of RX J0420 is low enough that it would not 
be remarkable in a radio pulsar. 

Second, the temperature of RX J0420 is the lowest 
of the INSs. Again, this might make sense if it is 
roughly the same age as the rest of the INSs but started 
with the lowest magnetic field. It would then have 
been heated the least, and would come closest to the 
"pristine" cooling of a non-magnetic neutron star. In 
this context, it is interesting that RX J1856.5— 3754 is 
cooler and less magnetized than RX J0720.4— 3125 (1.5 
vs 2.5 x 10 13 G), while appearing the younger one by 
kinematic age (Kaplan, van K crkwiik. & And ersonll2007l : 
iTetzlaff et all 120101 1201 ID. It would be interesting to 
measure the kinematic age of RX J0420 to compare it 
with the rest of the population. 

To view the evidence for field decay in a different way, 
we show in Figure [3] the blackbody temperature ver- 
sus characterist ic age for pulsars and the INSs (see also 
iZhu et all 120 111) . It is quite clear that the INSs are sys- 
tematically a factor of 5-10 older in characteristic age for 
the same temperature. If instead one uses the kinematic 
age, however, one sees that the difference is much smaller 
(for the two sources for which kinematic ages are avail- 
able). In the context of a picture in which the fields of 
INSs decayed, the main difference with pulsars induced 
by the initially much stronger field is thus that it leads 
to rapid initial spin down and long present periods (and 
thus long characteristic ages); the current temperatures 
are not as strong ly affected. Indeed, in the models of 
iPons et ail (|2009T ) ■ the heat generated by field decay is 
lost fairly rapidly. 

Third, given both the low temperature and low mag- 
netic field, RX J0420 largely follows th e empir ical 
temperature-magnetic field correlation from lKvK09l As 
discussed there, the origin of this relation (evolution- 
ary vs. surface physics) or even its overall integrity in 
the face of new data are not clear. It does seem to 
form an upper limit to the possible magnetic field of 
an INS, and even the rotating radio transient (RRAT) 
J1819— 1458 (possibly somewhat younger than the INSs, 
and with higher E as well) seems to roughly agree (based 
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Fig. 3. — Blackbody temperature (measured at infinity) ver- 
sus characteristic age for the INSs (circles) and rotation-powered 
pulsars (squares and upper limits); for all objects the color indi- 
cates the dipole magnetic field according to the scale at the right. 
The pulsars included here are both those that are visible in the 
ROSAT All-Sky Survey and those with dipol e field s > 10 13 G 
that have X-ray measurements. See IZhu et all H2011I) for a sim- 
ilar plot, although there they emphasize the variation among radio 
pulsars as a function of magnetic field, while here we emphasize 
primarily t he distinction between pulsars an d INSs. The data are 
taken from Kaplan & van Kcrkwiik (2009b|), with the addition of 
PSR B1916+14 UZhu et al.120091 ). PSR J1734-3333 HOlausen et~aTl 
2010fl . an updated temperature for PSR J1718-3718 UZhu et al.l 
20Tlll . an upper limit for PSR B1845-19 (Zhu et al. 2011, in 
prep), PSR J0726— 2612 (Kaplan et al., in prep.), an updated tim- 
ing solution for RX J0806.4— 4123 (van Kerkwijk et al., in prep.), 
and data from this paper. We only show data for sources where 
the blackbody can plausibly represent the majority of the surface 
(i.e., -Rbb > 3km, allowing for distance uncertainties). This is 
not meant to be a cooling curve, but instead to demonstrate the 
inconsistency in characteristic ages for the INS relative to the pul- 
sars at the same kT. The horizontal lines extending to the left of 
RX J1856.5-3754 an d RX J0720.4-3125 show the ranges of p lau- 
sible kinematic ages (Kaplan et al. 2007; Tctzlaff et al. 2011) for 
those objects, with the thick parts the mostly likely ranges. 



on lMcLaughlin et alll2007l ). 

Fourth, we did not confir m the tentative absorp- 
tion feature found by IH+04 although problems with 
the background subtraction meant we cannot refute it 
with confidence either. If it is true that RX J0420 
has no broad X-ray absorption feature, it would join 
RX J1856.5— 3754 (although as this object is often used 
for calibration assuming it emits like a blackbody, it is 
difficult to set confident limits). These are the two INSs 
with the lowest temperatures and the lowest magnetic 
fields, suggesting some relation between the presence of 
absorption features (or their energy) and either tem- 
perature or field strength (although a direct c orrelatio n 
of energy with field strength seems excluded; lKvK09). 
RX J0420 is also similar to RX J1856.5-3754 in its op- 
tical excess: both are reasonably well fit by Rayleigh- 
Jeans like pow erlaws, unlike the ot her INSs whose spec- 
tra are softer ( Kaplan et all 1201 ID . It is possible that 
the optical/UV spectral index is related to the mag- 
netic field, either directly through the magnetosphere 
(Tong ct al. 20 1l|) or indi rectly through shifting spectral 
lines ()Kaplan et ail 1201 ID : in that case the similarity of 
RX J0420 and RX J1856.5-3754 would be natural. 

Overall, our measurement firmly places RX J0420 as 
one of the INSs despite its short period, and moves us 
significantly closer to having a complete sampled of mea- 
sured spin-downs for that population. There are still a 
number of open questions to be answered via X-ray and 
multi-wavelength observations. Primary among these 
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is understanding the surface emission through consis- 
tent modeling of the spectra and lightcurves, and ideally 
with phase-resolved spectroscopy. Observations at opti- 
cal/UV wavelengths of the pulsed emission could make 
significant improvements in our understanding, by ty- 
ing the emitting areas at different wavelengths together 
and establishing the degree of surface inhomogeneiety. 
Finally, further kinematic ages would help greatly in 
constraining the coupled evolution of magnetic held and 
temperature. 



Based on observations obtained with XMM-Newton, 
an ESA science mission with instruments and contribu- 
tions directly funded by ESA Member States and NASA. 
DLK was partially supported by NASA through grant 
NNX08AX39G. Apart from the XMMSAS data reduc- 
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